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1. Introduction 

An atom moving across magnetic field is equivalent to an atom placed in perpendicular 
magnetic and electric fields. We consider the hydrogen atom moving in a magnetic field 
£?, strong enough to significantly squeeze the electron wavefunction. Quantitatively, 
the parameter 7 = 7iw c /(2Ryd) = £>/(2.35 x 10 9 G) is assumed large. Here, 
lu c = eB/m c c is the electron cyclotron frequency, and Ryd = m c e 4 /27i 2 is the ground- 
state energy of the field-free atom. 

Although only small values of 7 are available in laboratory, large values arc not 
uncommon in astrophysics. Spectra of some white-dwarf stars have been interpreted 
as produced by hydrogen at field strengths between 10 6 G and 10 9 G (Wunncr and 
Ruder 1987, Fassbindcr and Schwcizer 1996, and references therein). Neutron stars 
which are observed as radio pulsars reveal field strengths in excess of 2 x 10 8 G, and 
over half of them possess magnetic fields from 10 12 G to 2 x 10 13 G (Taylor et al 
1993). Absorption of radiation by strongly magnetized atomic hydrogen may have 
large effects on ultraviolet and X-ray spectra of the neutron stars, which are measured 
with modern space telescopes (Pavlov et al 1995). 

The physics of solid state presents another important field of application of 
quantum- mechanical calculations of strongly magnetized hydrogen atoms. Excitons 
and shallow impurities in semiconductors reveal hydrogen-like spectra with scaled 
values of lo c and Ryd. Such scaling offers a possibility to reach the regime 7 > 1 in an 
experiment (e.g., Elliott and Loudon 1960, Klaassen et al 1997). 

The non-moving hydrogen atom in a strong magnetic field was thoroughly studied 
in the past two decades (Ruder et al 1994, and references therein). Extensive tables 
of binding energies have been presented by Rosner et al (1984) and supplemented by 
Wintgen and Friedrich (1986), Ivanov (1988), Xi et al (1992) and Kravchenko et al 
(1996). Tables of oscillator strengths at various values of 7 have been published by 
Forster et al (1984); analytical fits to photoionization cross sections at 7 » 1 have 
been proposed by Potekhin and Pavlov (1993). 

The problem of a moving atom (or an atom in crossed fields) is much more 
complicated because of the absence of the axial symmetry. Much work has been 
done at low field strengths (e.g., Melezhik 1993), where the second-order perturbation 
theory was applicable (e.g., Braun and Solov'ev 1984, and references therein). The 
simplifying approximation of infinite proton mass, exploited in this regime, breaks 
down in strong fields because of collective motion effects studied in detail by Avron 
et al (1978), Baye and Vincke (1990) and Dippel et al (1994). In particular, so- 
called decentred states (with the electron localized mostly in the "magnetic well" 
aside from the Coulomb centre) are likely to be populated. These exotic states have 
been predicted by Burkova et al (1976) and studied by Ipatova et al (1984), Baye et al 
(1992), Dzyaloshinskii (1992) and Schmelcher (1993). As well as the usual "centred" 
states, the decentred states have an infinite discrete energy spectrum (Potekhin 1994, 
hereafter Paper I). Collective-motion effects on the centred states of the strongly 
magnetized hydrogen atom have been considered by Vincke and Baye (1988) and 
Pavlov and Meszaros (1993) in frames of the theory of perturbation. 

Completely non-pcrturbative results, covering both centred and decentred states 
as well as the transition region, were first presented by Vincke et al (1992) for binding 
energies and wavefunctions. In Paper I, additionally, oscillator strengths have been 
considered. Pavlov and Potekhin (1995, hereafter Paper II) studied spectral line 
shapes, and Potekhin and Pavlov (1997) calculated photoionization cross sections. 
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None of these numerical results, however, has been published in an easy-to-use form of 
tables or analytical expressions. This paper provides such expressions for the magnetic 
field strengths typical of neutron stars, 300 < 7 < 10 4 . This range is physically 
distinguished, since at weaker fields the transition region is strongly complicated by 
multiple narrow anticrossings (Vincke et al 1992). The relative simplicity of the 
spectrum at 7 > 300 facilitates analytical description. The upper bound, 7 <~ 10 4 , 
corresponds to the onset of non-negligible relativistic effects (Chen and Goldman 
1992). 

In the next section we recall the basic definitions and physical properties of a 
hydrogen atom arbitrarily moving in a strong magnetic field. In Section 3, we first 
present accurate analytical fits to binding energies, depending on the state of motion, 
for a number of bound states and various field strengths. Then we derive analytical 
approximations continuously depending on 7. As a by-product, simple and accurate 
approximations are obtained for binding energies of the non- moving atom at any 7 > 1. 
The obtained formulae are applied to evaluation of the maximum transverse velocity 
of the strongly magnetized atom. Section 4 is devoted to analytical approximations 
of quantum-mechanical sizes and main oscillator strengths of the atom. In section 5, 
an example of using the obtained expressions for calculation of absorption coefficients 
of strongly magnetized, hot hydrogen plasma is presented. 

2. Centred and decentred states: general description 

Motion of the hydrogen atom in a magnetic field can be conveniently described by the 
pseudomomentum (e.g., Johnson et al 1983) 

K = m p r p + m c r c - -B x (r e - r p ), (1) 

where the subscript i —e or i =p indicates electron or proton, respectively, 

Vi = ^[H tot , ri ] = -— V, - -^-A(n) (2) 

n mi m; L c 

is the velocity operator, rrii the mass, q p = —q c = e the charge, A(r) the vector 
potential of the field, and H to t the two-particle Hamiltonian operator. Gorkov 
and Dzyaloshinskii (1968) have shown that in the representation in which all three 
components of K have definite values the pseudoseparation of the centre-of-mass 
motion can be performed, that is, the relative motion can be described in terms of a 
one-particle Hamiltonian which depends on K. The expectation value of the velocity 
of the atom is Vjc-E, where E is the expectation value of the energy. 

It is convenient to describe the centred states of the atom using the relative 
coordinate = r c — r p as independent variable and the axial gauge of the vector 
potential, A(r) = \B x r. For the decentred states, the "shifted" representation 
(Gorkov and Dzyaloshinskii 1968) is more convenient. In the latter representation, 
the independent variable is = r — r p — r c and the gauge is A{r) = \B x (r — 
[(m p — m e )/mn] r c ). Here, r c = (c/eB 2 )B x K is the relative guiding centre and 
tor = mp + m c is the mass of the atom. 

Let us assume that B is directed along the z-axis. The z-component of the 
pseudomomentum corresponding to the motion along the field separates exactly from 
the Hamiltonian, giving the kinetic term Kz/2mn, while the transverse components 
K± produce non-trivial effects. Therefore we assume K z — and K± = K hereafter. 
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If there were no Coulomb attraction, then the electron Landau number n = 
0,1,2,... and the z-projection s of the angular momentum of the relative motion 
would be exact quantum numbers (since K is definite, the electron and proton do not 
possess definite z-projections of the angular momenta separately from each other - 
see Johnson et al 1983). In this case the transverse part of the wavcfunction could be 
described by a Landau function $„ s (r^), where is the projection of in the 
(xy)-plane and s is defined in the shifted reference frame (e.g., Paper I). The energy 
of the transverse excitation (with the zero-point and spin terms subtracted) is 

E ns = [n + {me/ nip) (n + s)]huj c . (3) 

A wavefunction ip K of an atomic state |k) can be expanded over the complete set 
of the Landau functions 

^\r^) = Y J ®ns{r { 2 ) )9^ K {zl (4) 

ns 

where rj — or 1 indicates the conventional or shifted representation, respectively (a 
generalization to continuous rj in Paper I proved to be less useful). The adiabatic 
approximation used in early works (Gorkov and Dzyaloshinskii 1968, Burkova et al 
1976) corresponds to retaining only one term in this expansion. 

A bound state can be numbered as \k) = \nsvK), where n and s relate to the 
leading term of the expansion (4) , and v enumerates longitudinal energy levels 

E\ SV {K) = E K -E^ S (5) 

and controls the z-parity: ($} S - K {— z) = (— 1)" g n v l- K {z). This way of numbering is 
conventional for the non-moving atom at 7 > 1. The states v = are tightly bound 
in the Coulomb well, while the states v > 1 are hydrogen-like, with binding energies 
below 1 Ryd. For a moving atom, this way of numbering remains unambiguous at 
7 > 300, in spite of the fact that there may not exist an obvious leading term of (4) 
in this case (Paper I). 

The inequality E K < determines truly bound states, as opposed to the ones 
subject to autoionization. In particular, all states with n ^ belong to continuum at 
7 > 0.2 and will not be considered hereafter. 

Since the transverse factors & ns in (4) arc known analytically, only the one- 
dimensional longitudinal functions are to be found numerically. An algorithm 
which is most efficient at 7 » 1 has been described in Paper I. At small 
pseudomomenta K, the states v = remain tightly bound and centred, the average 
electron-proton displacement x being considerably smaller than r c . For the hydrogen- 
like states v > 1, however, x is close to r c at any K. 

According to the second-order perturbation approximation at small K, the 
absolute expectation value of the velocity v = 8E K /dK in a bound state \k) equals 
K/M^sv, where M^ sv is the effective "transverse" mass (Vincke and Baye 1988, Pavlov 
and Meszaros 1993). M^ sv always exceeds tor, and it is the greater the stronger the 
field and the higher the considered atomic level. 

The larger K, the greater is the distortion of the wavefunction towards r c , 
caused by the motion-induced electric field in the co-moving reference frame. The 
perturbation approximation becomes increasingly inaccurate, until near some critical 
value K c a transition to the decentred state occurs, and the character of the motion 
totally changes. With further increasing K, the transverse velocity decreases and 
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tends to zero, while the electron-proton separation increases and tends to r c . Thus, for 
the decentred states, the pseudomomentum characterizes electron-proton separation 
rather than velocity. 

The shifted (77 = 1) adiabatic approximation becomes fairly good at K » K c . 
At very large K the longitudinal functions become oscillator-like, corresponding to a 
wide, shallow parabolic potential well of a depth ~ e 2 /r c (Burkova et al 1976). For a 
fixed v, this limit is reached at K ^> [y + ^) 2 h/aB, where oq is the Bohr radius. Still 
at arbitrarily large K there remain infinite number of bound states with high values 
of v whose longitudinal wavefunctions are goverened by the Coulomb tail rather than 
by the parabolic core of the effective one-dimentional potential (Paper I) . 

The decentred states of the atom at K > K c ~ 10 2 au have relatively low binding 
energies and large quantum-mechanical sizes, I ~ K/j au; therefore they are expected 
to be destroyed by collisions with surrounding particles in the laboratory and in 
the white-dwarf atmospheres. In neutron-star atmospheres at 7 > 10 3 , however, the 
decentred states may be significantly populated (Paper II). This necessitates inclusion 
of the entire range of K below and above K c in the consideration. 



3. Binding energies 

3.1. Dependence of the energies on the pseudomomentum at selected field strengths 

We have calculated binding energies of the hydrogen atom moving across the strong 
magnetic field at 7 = 300, 600, 1000, 2000, 3000 and 10 000 for several lowest tightly- 
bound and hydrogen-like states, using the technique described in Paper I. At each 
value of 7 and for each state, the calculations have been performed at K — and at 
about 50-100 values of K from K < 10 au to K > 10 4 au, approximately equidistant 
in log if but with additional points near avoided crossings. The calculated energies 
have accuracy of 3-5 digits. 

In applications, however, one usually has to deal with a distribution of atoms over 
more or less broad band of values of the pseudomomentum K, and to calculate the 
observable quantities by averaging over K. This makes it highly desirable to have an 
analytical approximation of the if-dependence of the energies, E(K). Lai and Salpeter 
(1995) were the first to present an analytical fit to E(K), which was rather accurate 
for the ground state at K < K c but could not be applied to excited or decentred 
states. 

We describe the longitudinal energy (5) by the formula 

, F \\ ( K \\ = E ( nX(K) E { n} u (K) 

I nsvV )\ l + ( K / Kc y/a "•" l + (Jf c /A-)l/a' [ > 

The two-term structure of (6) is dictated by the necessity to describe the two physically 
distinct regions of K below and above K c . The parameter a has the meaning of the 
width of the transition region near K c in logarithmic scale of pscudomomcnta. As 
noted in Paper I, for the tightly-bound states K c is close to (2muEnJ„) 1 ^ 2 , where 
En)v = —Ensi>(0). We write K c = q^m^En^) 1 ^ 2 and treat go as a fitting parameter. 

Intricate structure of the region of avoided crossings (see Vincke et al 1992) 
complicates its accurate analytical description. We have chosen to keep our formulae 
simple at cost of decreasing accuracy near these crossings. 
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For the tightly-bound states, we parametrize the functions E^ (K) as follows: 

E oso( K ) = E oll - ~777oT' ( 7 ) 



Eg> (K) = 2[rj+ri /z + qin\ ' Ryd, (8) 

where r* = r c /a-Q = K/(j au), q\ and q 2 are dimensionless fitting parameters, and 
m ff is the effective mass which is close to (but not necessarily coincident with) the 
transverse effective mass M^ sv obtained by the perturbation technique. 
In the considered range of 7, the parameter qi can be approximated as 

= f lg( 7 /300) if s = 0, 
1 [ 0.5 otherwise. 

Optimal values of the other parameters are listed in table 1. The last column presents 
the root-mean-square (rms) difference ue between the computed and fitted energies. 
Maxmum errors of the fit ( < 3a e) occur near the avoided crossings. 

Binding energies of the hydrogen-like states are approximated by the same formula 
(6) but with slightly different expressions for E^> and E^ 2 \ For these states, M^ sv 
exceeds toh by orders of magnitude, and the perturbation method fails already at 
small values of K (Pavlov and Meszaros 1993), which renders the notion of transverse 
mass practically useless for the fitting. Thus we consider m e ff as effectively infinite 
and put EqU = E^l (y > 1). Furthermore, the transition region is not well defined, 
and therefore K c and a lose their clear meaning and become mere fitting parameters. 

The function E^ (K) that describes the longitudinal energy at large K is now 

E£l(K) = I (2 Ryd)- 1 [rl + {2v + l)rl /2 + q 2 r*] ^ + (9) 

where r* and E^ have the same meaning as before. The first and second terms in 
the square brackets ensure the correct asymptotic behaviour (Paper I). In this case, 



,2 , 3/2 



-1/2 



<Z2 



V 2 - 1 (odd V) 

v 2 + 2"/ 2 lg(7/300) (evenj/). 



Optimal values of the parameters go and a are listed in table 2. As well as in table 1, 
the last column presents rms errors which are several times smaller than the maximum 
errors near anticrossings. 

In both tables 1 and 2, only truly bound (not autoionizing) states are considered. 
For example, all states with s > 0, v > belong to continuum at 7 > 673, therefore 
table 2 does not contain entries for them at 7 = 1000 and higher. 



3.2. Two-dimensional approximations 

Equations (6)-(9) help us to derive approximations of the binding energies as functions 
of two continuous arguments 7 and K. For this purpose, we replace the numerical 
parameters listed in tables 1 and 2 by analytical functions of 7. 

One of these parameters — the longitudinal energy of the atom at rest E^ - 
has an independent significance. For this reason, we have constructed an accurate fit 
to E^ in a possibly widest range of 7 values. For the tightly-bound states, we have 

4° ( 7 )/Ryd = In (cxp [(1 + s)- 2 ] + Pl [ln(l +P2Vl)} 2 ) +P3 [ln(l +P47 P5 )] 2 • (10) 
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Table 1. Parameters of the analytical approximation (6)— (8) for the energies of 
tightly-bound states |0s0) 



s 


7 


(Ryd) 


lg(m e ff/m. H ) 


90 


a 


92 


a E (Ryd) 





300 


10.722 


0.009 


0.859 


0.001 


0.102 


0.028 




600 


13.210 


0.042 


0.811 


0.107 


0.157 


0.040 




1000 


15.325 


0.072 


0.823 


0.117 


0.189 


0.025 




2000 


18.610 


0.141 


0.850 


0.178 


0.233 


0.018 




3000 


20.770 


0.175 


0.873 


0.191 


0.244 


0.017 




10000 


28.286 


0.319 


1.019 


0.173 


0.275 


0.027 


1 


300 


7.669 


0.161 


0.963 


0.132 


0.115 


0.026 




600 


9.607 


0.269 


1.060 


0.093 


0.160 


0.021 




1000 


11.277 


0.369 


1.147 


0.060 


0.176 


0.024 




2000 


13.904 


0.578 


1.195 


0.122 


0.215 


0.016 




3000 


15.649 


0.701 


1.202 


0.147 


0.235 


0.014 




1UUUU 


91 8^0 


n QAA 


1 ^37 


fl 9Q& 

u.zyo 






2 


300 


6.450 


0.304 


1.184 


0.030 


0.120 


0.017 




600 


8.142 


0.497 


1.197 


0.081 


0.181 


0.014 




1000 


9.610 


0.643 


1.262 


0.074 


0.195 


0.014 




9nnn 
zuuu 


1 1 0.^7 

ii.yo i 


u.yo± 


1 9Q1 


n 1 97 

U. 1Z I 


o 9^n 

U.ZOU 






OUUU 




i.uyo 


1 ^90 


n i w 

U. 10o 


o 9zin 


fl 099 

u.uzz 


3 


300 


5.734 


0.466 


1.263 


0.039 


0.122 


0.015 




uuu 


7 974 


n 701 


1 97^ 
± .z / o 


n ns9 


u. loo 


n ni 9 

U.U1Z 




1UUU 




fl £Q7 

u.oy ( 


1 ^A7 


fl 00.0 

u.uyu 


90A 


U.Ulo 




2000 


10.755 


1.252 


1.403 


0.131 


0.232 


0.019 




3000 


12.191 


1.451 


1.457 


0.154 


0.240 


0.026 


A 
1 


oUU 


^ 9/1 Q 


U.OlO 


l.ooU 


U.UOU 


n i 9Q 

U. IZo 


U.Ulo 




600 


6.676 


0.892 


1.342 


0.095 


0.194 


0.011 




1000 


7.929 


1.124 


1.437 


0.096 


0.211 


0.017 




2000 


9.933 


1.555 


1.544 


0.114 


0.229 


0.016 


5 


300 


4.877 


0.755 


1.391 


0.058 


0.128 


0.012 




600 


6.227 


1.086 


1.393 


0.107 


0.199 


0.012 




1000 


7.413 


1.354 


1.545 


0.130 


0.229 


0.010 


6 


300 


4.589 


0.888 


1.448 


0.062 


0.123 


0.013 




600 


5.874 


1.281 


1.441 


0.121 


0.207 


0.013 




1000 


7.004 


1.668 


1.587 


0.107 


0.210 


0.018 


7 


300 


4.355 


1.021 


1.504 


0.070 


0.132 


0.013 




600 


5.585 


1.480 


1.473 


0.139 


0.213 


0.014 



The parameters pi — p$ depend on s; they are presented in table 3. This fit is accurate 
to within 0.1-1% at 7 = 10 _1 — 10 4 , and it also provides the correct limits at 7 — ► 0. 

For the hydrogen-like states, we use the asymptotic result (Haines and Roberts 
1969) 

« = 7^. where { " = ^ V" ^ "\ (H) 

(n + (5) 2 ( n — v/2, ~ (In 7) (even z/). v y 

We have obtained the following fits to the quantum defect <5: for odd 1/, 



<5( 7 ) = (a, y + fe,V7 + 0.077 7 )- 1 , 



(12) 



■55 



Table 2. Parameters of the analytical approximation (6), (9) for the energies of 
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300 


0.08920 


2.435 


0.8688 


0.00054 






600 


0.09068 


4.328 


0.7156 


0.00016 






1000 


0.09167 


5.237 


0.7205 


0.00018 






2000 


0.09294 


7.419 


0.6593 


0.00010 






3000 


0.09362 


8.825 


0.6237 


0.00016 






10000 


0.09542 


13.43 


0.5906 


0.00026 


1 


1 


300 


0.9407 


2.109 


0.6794 


0.0010 






600 


0.9640 


3.553 


0.7029 


0.0024 


1 


2 


300 


0.5138 


1.930 


0.6417 


0.0038 


2 


1 


300 


0.9223 


2.421 


0.6553 
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Table 3. Parameters of the analytical approximation (10) for the energies of tightly- 



bound states 


|0s0) at 


10 1 < 


y < 10 . 




s 


Pi 


P2 


P3 


P4 


P5 





15.55 


0.378 


2.727 


0.3034 


0.4380 


1 


0.5332 


2.100 


3.277 


0.3092 


0.3784 


2 


0.1707 


4.150 


3.838 


0.2945 


0.3472 


3 


0.07924 


6.110 


4.906 


0.2748 


0.3157 


4 


0.04696 


7.640 


5.787 


0.2579 


0.2977 


5 


0.03075 


8.642 


6.669 


0.2431 


0.2843 


6 


0.02142 


9.286 


7.421 


0.2312 


0.2750 


7 


0.01589 


9.376 


8.087 


0.2209 


0.2682 



where a v s=a 1 and b v ~ 2; and for even v, 



n -1 



(13) 



| and b v ~ |. Accurate values of a v and b v are given in table 4. At 



where a v 

1 < 7 < 10 4 , rms errors of (12) lie within 3 x 1CT 4 , and those of (13) within 10" 



Table 4. Parameters of the analytical approximations (11)— (13) for the energies of 
hydrogen-like states |00i^) at 1 < 7 < 10 4 . 



V 


1 


2 




4 




6 


a v 
K 


0.785 
1.724 


0.578 
0.765 


0.901 
1.847 


0.631 
0.717 


0.970 
1.866 


0.660 
0.693 



The parameters m c ff, a and qo in (6)-(9) that determine if -dependences of the 
energies can also be replaced by analytical functions of 7. Let us start with the 
tightly-bound states (y = 0). For the effective mass, we have 

m eff (7) = m H [1 + (7/7o) Co ] , (14) 

where the power index c and the value 70 (roughly corresponding to the onset of 
strong coupling between internal and ccntre-of-mass motions of the centred atom) 
depend on the quantum number s and are given by 



c = 0.937 + 0.038s 158 and 7o = 6150 (l + 0.0389s 3/2 ) 1 + 7.87s 3/2 
For the critical pseudomomentum, we write 

q a = K c /y/2m H E(°) = d + ln(l + 7 / 7l ). (15) 

The parameters ci and 71 take on the values c\ — 0.81,1.09,1.18,1.24 and 71 = 
(8.0,3.25,2.22,1.25) x 10 4 for s = 0,1,2,3, respectively For s > 4, we put 
ci = 0.93 + 0.08s and 71 = 6500. The remaining parameters can be replaced by 
simple expressions, a = 0.053 ln(7/T50) and q 2 = 0.158 [ln((l + 0.1s) 7 /215)] 2 / 5 . 

Now let us turn to the hydrogen-like states. For odd states, we have, 
approximately, q a — (^ 5 / 4 7/170) 9 and a — 0.66 + ^/20, whereas for even hydrogen- 
like states q = v ^/7/1200 and a = 0.66. 




Figure 1. Energy spectrum of the hydrogen atom moving across a strong magnetic 
field. Upper panels: tightly-bound states (v = 0); lower panels: hydrogen-like states 
of the manifold s = 0. Numerical values (0) are compared with the analytical 

approximations of section 3.2 ( ) and of Lai and Salpetcr (1995) ( , upper 

panels). Triangles (A) mark the limit of validity of the perturbation formalism 
according to equation (3.8) of Lai and Salpeter; large circles (O ) point the present 
analytical approximation for the critical pseudomomcntum K c . 



These approximations are not so accurate as those provided by tables 1 and 2, but 
their advantage is that they may be used at arbitrary 7 in the range considered. In 
figure 1 they are compared with our numerical results and with the fitting formulae 
of Lai and Salpeter (1995). The figure demonstrates that the present approximations 
are valid at any K from to infinity. Noticeable discrepancies between our fitted and 
calculated data occur only in narrow ranges of K near anticrossings. 
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3.3. The largest transverse velocity 

As an example of application of the above formulae for binding energies, let us estimate 
the maximum velocity of the atom, t> max = max \dE/dK\. It can be alternatively 
interpreted as the maximum transverse electric field £ mov — v mllx B/c that could be 
applied to an atom at rest. A stronger electric field £ > £ m0 v forces the atoms 
to move with velocities around the drift velocity of free charges in crossed fields, 
Adrift = c£ x B/B 2 , provided that £ < £ c = B = 137.0367 au ( m conventional 
units, £ c = 2.998 x 10 4 [B/G] V nT 1 ). Still higher electric field, £ > £ c , cannot be 
counterbalanced by motion, hence it causes Stark ionization. 

A numerical evaluation of v max requires multiple calculation of derivatives dE K /dK 
in the most complicated region K <~ K c . Thus the analytical approximations can be 
most helpful here. A reasonable approximation is simply w max .appr = d\E\/dK at 
K = K Cl where E is given by (6) and K c by (15). Table 5 presents t> m ax,appr obtained 
using this approximation for the lowest tightly-bound states along with v max ,LS given 
by equation (3.30) of Lai and Salpeter (1995) and with w max evaluated numerically. 
The values listed in table 5 (in atomic units of velocity, 1 au= 2188 km s _1 ) can be 
translated into those of the critical electric field, £ mov — ■yv max [au] ■ (5.14 x 10 11 V 
m- 1 ). 

Table 5. The largest transverse velocities (in au) in the lowest states OsO): numerical 
values t>max compared with the present analytical approximation t> ma x,appr and with 
the approximation t) max ls °f L a i an d Salpeter (1995). 



s 
7 




300 




1000 




3000 




10 000 


1 

3000 


2 

3000 


3 

3000 


%ax 

^max.appr 

fmax.LS 


0.0588 
0.0622 
0.0812 


0.0479 
0.0467 
0.0551 


0.0372 
0.0367 
0.0374 


0.0253 
0.0279 
0.0228 


0.0240 
0.0232 
0.0174 


0.0198 
0.0187 
0.0124 


0.0173 
0.0154 
0.0100 



4. Geometrical characteristics and radiative transitions 

4-1. Atomic sizes and dipole moments 

Geometrical characteristics of an atom play an important role in distribution of 
atoms over quantum states in a plasma and in their contribution to the plasma 
absorption coefficients, since a "size" of an atom may be used to approximately 
evaluate effects of destruction of atoms caused by random charge distribution in the 
plasma (e.g., Potekhin 1996). The if -dependence of rms sizes is complicated and can 
be non-monotonous near anticrossings. However, the sizes usually need not be known 
with high precision, that relieves the problem of fitting. The accuracy level of the 
approximations presented in this section is typically several percent. 

At K = 0, the atom is axially symmetric, and its rms sizes along the Cartesian 
coordinates can be approximated as I x q = l y o ~ ae \/(s + and 

k„« {l/v / 2 + l/ln[ 7 /(l + .s)]}(Ryd/i : ;( )) 1 / 2 a B (v = 0) (16) 

U«(1.6Ryd/^°))aB {v > !)• (17) 

Let us consider an atom moving along Oy. Both transverse sizes of the electron 
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"cloud" remain approximately independent of K: l x w l y rj l x0 . However, the atom 
acquires a constant dipole moment d — e(r p — r e ) proportional to the mean proton- 
electron separation x = \ (r e — r p )\. This separation is always smaller than r c , and it 
approaches r c at K 3> K c . With inaccuracy up to 10%, at 7 > 300, 

- l + 0.015 7 2 \/TT^(£ (0) r 4 r 1 [l + {K/K c f' c 
- 1 + 0.004 7 2 (E^ (K))- 4 ] ~* [l + (K c /K)^ a ] ~* , (18) 



x/r c 



where E<- 2 \ K c and a are defined above. 

The size of the electron cloud along the field is also affected by the motion. It can 
be described by the formulae 



h — ho 



[l + {l-m K /m eS )(K/K c y 



,1/2 



+ 



lz2 



1 + (Kc/K) 1 ^ 



(!/>!). 



(19) 
(20) 



1 + (K/Kc) 1 /* 

lz = (lla + II2) 
Here, l z0 is the value at K = given by (16), (17), and 

l z2 = ^7TTj2 [ r 3 + ( 4 .3 + 7 Z ,> 2 ] 1/4 
has the correct asymptotics at K — > 00 (Paper I) . 

1 /2 

In figure 2 the average size of the atom, l(K) = \x 2 + l x + l y + I 2 ] , expressed 
through the above formulae, is compared with values calculated numerically. On the 
left panel (at 7 = 600), the strong deviations of the numerical values from the fit 
for the states 1 003) and 1 1 1) at if ~ 10 2 au are caused by their anticrossing, which 
occurs shortly before the level (011) enters continuum. At 7 = 2000 (right panel), this 
level belongs to the continuum at arbitrarily small K, so it does not (anti)cross truly 
bound levels. 



4-2. Oscillator strengths 

In this section we consider those oscillator strengths / that dominate photoabsorption 
of polarized radiation by ground-state hydrogen atoms at large 7. The polarization is 
assumed circular (right, for which we will use superscript '+', or left, '— ') or linear, 
longitudinal ('||') with respect to the static magnetic field. 

At K = 0, the left-polarized radiation cannot excite the ground-state atom 
(/~ = 0), while right and longitudinally polarized radiation is absorbed mainly via 
transitions to the states |010) and |001), respectively. The corresponding oscillator 
strengths have been computed and tabulated by Forster et al (1984). With an accuracy 
of 1-2%, at < 7 < 10 4 , they are reproduced by a single five-parameter formula 

&,(0)=fl- °- 584 ) 1 + M37 - (21) 

For /oio(0)j the parameteres m take on the values u\ = 12, 112 = 1.43, 113 = 9.8 x 10~ 5 , 
u A = 1.585 and u 5 = 0.713. For /J 01 (0), we have m = 2.64, u 2 = 1.076, m 3 =6x 10~ 6 , 
u 4 = 0.247 and u 5 = 0.381. 

For the moving atom, in the restricted range 300 < 7 < 10 , we put 

f+ (K) - f+ (0) 1 ~ < K l K <y + 2 ( m */ m p) (22 ) 



GO 




|000) 



i i i i i i i i i i i i i i i i i i i i i i i i i i 

1 2 3 1 2 3 

lg K (a.u.) lg K (a.u.) 

Figure 2. Comparison of calculated geometrical sizes of the atom (0) with the 
analytical approximations of section 4.1 ( ). 



with a = 1.28 - 0.267 ln(l + 7/240), 6=1 + 3/[l + ln 2 (l + 7/90)], a' = 0.012 [1 + 
In 2 (1 + 7 /90)], and 



flM = fU0) eM-(a'K/K c f] + ^ f? 1 *^ , (23) 

1 + U.O^K c /K 

with a' = 0.8771n(13100/7), V = 0.89 - 7 /17000 and (3 = 0.61 (1 + 2410/7) 3 / 2 . The 
second parts of (22), (23) ensure the correct large- if limits (2m e /m p and 1 respectively, 
cf. Paper I). 

The radiative transitions forbidden for the atom at rest because of the 
conservation of the angular-momentum projection become allowed for the moving 
atom. In particular, the moving ground-state atom can absorb left-polarized radiaton. 
Oscillator strengths of such transitions are significant only at K of the order of 
K c ~ 10 2 au. Therefore we derive for them fitting formulae accurate to <~ 10% 
in this range of K and do not attempt to fit the complicated behaviour they show 
outside this range, where they are orders of magnitude smaller (Paper I). 

The transition to a state |0s0) presents the dominating absorption channel for 

circular polarization in a spectral range — E^°l_ 10 < Hlu < — Eq°J (where 
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hoj is a photon energy). For the right polarization, we put 

+ _ 0.012 {K/K c r(l - K/K c ) 
Ioso[K ] ~ 1 + 11 ln(l + ( 7 /3300)2) at A < K Cl s > 2, (24) 

and zero at K > K c ; for the left polarization, 

hso{*)- 2 >[l + (K/K c )*l>+»] [ ~ h { ' 

Although approximations (22)-(25) are rather crude, particularly owing to the 
anticrossings, their accuracy may be still sufficient for astrophysical applications, as 
will be demonstrated in the next section. 



5. Spectral line shapes 

As an application of the above fitting formulae, let us consider bound-bound 
absorption spectrum of hydrogen under the conditions typical for neutron star 
atmospheres (Pavlov et al 1995): density p > 10~ 2 g cm~ 3 , temperature T <~ 
10 5 — 10 6 K, and magnetic field strength B ~ 10 12 — 10 13 G. Such absorption spectra 
have been studied in Paper II. Neglecting the Dopplcr and collisional broadening but 
taking into account the most important magnetic broadening, one obtains an average 
partial cross section of an atom with respect to absorption of polarized radiation with 
frequency u> via some specific transition \k) — > \k'} as the sum cr(u>) = Xi^O^) over 
the roots K^lu) of the equation E'(K) - E(K) = Huj, where E'(K) (E(K)) is the 
energy of the final (initial) state of the atom, and 

n 47r 3 e 2 

<n{w) = C' 1 K^dKjdivlw'iK^expi-EiK^/kBT}/^). (26) 

JTIqC 

Here, f(K) is the oscillator strength for the given transition and polarization, ks is 
the Boltzmann constant, 

poo 

C w = 2tt Kw(K)exp[-E(K)/k B T]dK (27) 
Jo 

is a normalization constant, and w'{K) (w{K)) is the occupation probability of the 
final (initial) atomic state in a plasma with a number density of electrons n e . These 
occupation probabilities can be estimated as w ~ exp[— (47r/3)n (47) 3 ], where I is 
the rms radius of the atom. This expression is pertinent to calculation of bound- 
bound absorption considered here. For calculation of thermodynamic properties of the 
plasma, however, one should use, instead of these "optical" occupation probabilities, 
"thermodynamical" ones, which are generally larger (Potckhin 1996). 

In order to take into account induced emission, cn{uo) should be multiplied by 
(1 — e -Hu/ k BTy co nisioiial broadening can be taken into account by convolution 
of (Tj(w) with the Lorentzian profile characterized by the width T(Ki(uj)) (Paper II). 
Since this type of broadening has only marginal significance compared to the magnetic 
broadening in the neutron star atmospheres, we employ a simple order-of-magnitudc 
estimate: 



T(K) w r n o a 3 (fc B T/Ryd) 1/6 (l + 2r^ /6 ), 



(28) 
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where T w 15 au for transitions to the state |001) (that mainly determine absorption 
of radiation polarized longitudinally) and To w (68/7) au for transitions to the state 
1 010) (responsible for the main absorption peak of circular polarization). 

Some of the typical absorption profiles obtained numerically in Paper II are 
represented by full lines in figure 3. The spectral range is shown that is relevant 
to interpretation of spectral observations of neutron stars with the X-ray telescope 
on-board ROSAT satellite (e.g., Pavlov et al 1995). Approximate profiles, obtained 
using the analytical approximations of section 3.2 for E(K), the formulae of section 4 
for l(K) and f(K), and the estimate (28) for T(K), are shown in figure 3 by broken 
lines. One can see that they correctly reproduce gross features of the spectral shapes, 
and thus the proposed approximations are suitable for using them in theoretical models 
of neutron-star atmospheres compatible with contemporary observational data. The 
figure demonstrates also the importance of the decentred states: for example, the 
absorption of longitudinally polarized radiation at hu; < 100 eV is produced solely by 
the states of motion with K > K c . 

In reality, observed spectra arc influenced not only by the bound-bound 
photoabsorption shown in figure 3 but also by bound- free transitions. For moving 
atoms, bound-free cross sections were calculated by Potekhin and Pavlov (1997). 
Under conditions typical of atmospheres of cooling neutron stars, maximum bound-free 
photoabsorption turns out to be of the same order of magnitude as the bound-bound 
one, but it is shifted to higher photon energies. In the spectral range shown in figure 
3, the bound-free transitions are generally less important. 



6. Conclusions 

We have obtained analytical approximations of binding energies, geometrical sizes 
and main oscillator strengths of radiative transitions of the hydrogen atom moving 
across a strong magnetic field. These approximations can be also applied to the 
hydrogen atom in crossed electric and magnetic fields, since the latter problem reduces 
to the former one with the effective pseudomomentum K = m\ic£ /£ c , or equivalcntly 
K [au] = 8.4 (£/V m" 1 ) (B/G)" 1 . 

Binding energies are the most important quantities in many applications, and 
for that reason we have presented not only fitting formulae analytically depending 
on 7 and K (section 3.2), but also considerably more accurate if-dependences at 6 
selected values of 7 (section 3.1). Atomic sizes (section 4.1) play important role in 
distribution of atoms over quantum states in a plasma and in their contribution to 
the plasma absorption coefficients. For example, a size of an atom may be used to 
evaluate effects of "unbinding" of electrons caused by random charge distribution in 
the plasma. For non-magnetized hydrogen plasmas, an approximate treatment of these 
effects was revised recently (Potekhin 1996); for the strong magnetic fields analogous 
work is under way. The approximations of oscillator strengths (section 4.2), along 
with those of the energies and sizes, facilitate calculations of absorption spectra of 
strongly magnetized, partially ionized hydrogen plasmas. Eventually, the analytical 
estimates of 7- and if -dependences of the binding energies, atomic sizes and transition 
rates will help to generalize previously developed models of fully ionized atmospheres 
of magnetic neutron stars (Shibanov et al 1992) to the more realistic case of partially 
ionized atmospheres. 
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Figure 3. Comparison of numerically calculated ( ) and approximate ( ) 

spectral shapes for bound-bound absorption of right (+), left (— ) and longitudinally 
(||) polarized radiation by hydrogen atoms in thermal equilibrium with plasma under 
the conditions typical for atmospheres of magnetic neutron stars: k^T = 1 Ryd 
(T = 1.58 x 10 5 K), p = 0.01 g cm" 3 , and 7 = 10 3 (B = 2.35 x 10 12 G). 
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